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We develop a theory for micro- Raman scattering by single and coupled two-donor states in silicon. 
We find the Raman spectra to have significant dependence on the donor exchange splitting and the 
relative spatial positions of the two donor sites. In particular, we establish a strong correlation 
between the temperature dependence of the Raman peak intensity and the inter-donor exchange 
coupling. Micro-Raman scattering can therefore potentially become a powerful tool to measure 
interqubit coupling in the development of a Si quantum computer architecture. 
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Quantum computing has attracted wide-spread inter- 
est in recent years Among many proposed schemes 
for quantum computer (QC) hardware, a particularly 
promising architecture is the nuclear spin based Si QC 
based on donor arrays 2], which has prompted a great 
deal of activities to study both theoretically and exper- 
imentally the fabrication, control, and measurement of 
such Si-based devices |llillli,0lEllE3,E3- Currently 
the fabrication of donor arrays in Si has been pursued 
using two complementary approaches. The "bottom-up" 
approach uses low temperature MBE and STM to build 
the array atom by atom 0, while the "top-down" ap- 
proach uses ion- implantation |l2l |. In both approaches, 
it is crucial to know the precise positions of the donors 
buried inside the host Si lattice, and to measure the ex- 
change coupling between neighboring donors. 

The most important (and daunting from the fabrica- 
tion perspective of building a Si QC) physical quantity 
underlying the Si QC architecture is the exchange entan- 
glement of the neighboring donor sites - this exchange 
coupling must be large enough for two-qubit operations 
to be possible [1,0,SM0- Right now, however, there is 
no proposed method for directly estimating the exchange 
entanglement in Si, and the experimental fabrication of 
Si QC architecture is essentially based on the hope of 
appreciable entanglement in the fabricated devices. In 
this Letter we not only propose a novel and promising 
quantitative method for obtaining the exchange coupling 
in Si donor states based on the easily available technique 
of electronic micro- Raman spectroscopy 13, 14], but also 
demonstrate through quantitative calculations the feasi- 
bility of our proposed technique as a powerful diagnostic 
tool for studying buried donor positioning and exchange 
coupling measurement. Electronic Raman spectroscopy 
has been widely used to study electronic transitions for 
which one-photon absorption |1ItI | is forbidden. Further- 
more, micro- Raman spectroscopy, with the laser spot size 
in the order of /xm, has the required spatial resolution to 



provide valuable information on the environment of single 
donors, pair orientation, and the strength of donor-pair 
coupling and electronic entanglement. In the following, 
we develop a theory for micro-Raman scattering from 
donor electrons in Si, focusing on the situation where 
photon polarization and sample temperature can be pre- 
cisely tuned. We show that the electronic micro-Raman 
spectra have important signatures of donor locations and 
exchange coupling. We also discuss specific experimental 
conditions when single donor (or single pair) detection 
might be possible with this technique. 

In Raman scattering (i.e. inelastic light scattering) 
experiments, an incident laser photon of wave vector k^, 
frequency lol and polarization ff^ scatters into a Raman 
photon with k^, ujr, ffn, while the system with which the 
radiation interacts undergoes a transition from an initial 
state |0) of energy Eq to a final state |/) of energy Ef. 
The differential Stokes- Raman scattering cross section is 
given by 
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where ujfo ~ [Ef — Eo)/h and to — cjl ^'-^r. is the Raman 
shift. The matrix element M is written as 
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where Eq is the direct band gap at the conduction 
band minimum and p is the momentum operator. In 
Si, the conduction band has six degenerate minima. The 
lowest donor-electron energy states are a ground-state 
singlet lS{Ai), and excited triplet lS'(r2) and doublet 
^S{E) components, consistent with the Td site symme- 
try of the substitutional donor [13. The notation lS{j) 
refers to a j-symmetry state obtained from IS-like hy- 
drogenic envelopes, which can be written as |^ '0i(r) = 
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Em=i"m^^mW<?^'mW: where (/-^(r) 
Bloch wavef unctions, J)j(r) are the hydrogenic envelopes, 

aji characterize the point group representation, and /i 
runs over the six conduction band minima k^. Taking 
this muitivaiiev feature into consideration, Eq. |^ can 
be simplified 
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where TZ is the resonance enhancement factor which, for 
h{ujL — w_r) <C Eg, can be expressed as 7^ = Eq/[Eq — 
(Tiwl)^], my and m_L are the Si conduction band eS^ective 
masses in units of the free electron mass m, and fj^ is the 
unit vector along the axis of the /ith valley. 

Electronic Raman scattering in Si:P is dominated by 
the valley-orbit transition of \S{Ai) iS{E) states for 
single donors (the lS{Ai) 1S{T2) transition is not 
Raman-active) 0|. The shift and width of the single 
donor Raman peak can provide donor environment in- 
formation, since the donor states are sensitive to the 
local strain fields, the local electric fields from ionized 
acceptor-donor pairs or surface charge for donors near 
the surface, and interactions with other neutral donors 
and acceptors. 

To have single donor sensitivity, a micro- Raman exper- 
iment needs sufficiently large Raman cross section, which 
is possible through the resonant Raman effect. The Ra- 
man cross section can be estimated neglecting the exci- 
tonic effects . In this case the divergence in the reso- 
nance enhancement factor TZ is limited by the energy level 
dispersion in the valence band, since the valence bands 
are not extremal at the conduction band minima. As a 
result, TZ can take on a maximum value of TZ^ « 1000, 
leading to a Raman cross section da/dVtji « 2 x 10^^ 
cm^/sr. For example, for a laser power of 1 mW the 
photon flux is 1.6 X 10^^ photons/s • cm^. This leads to 
a count rate of about 35 s~^ per donor, small but work- 
able. In addition, at T = lOK and for the incident laser 
power of 1 mW in a 1 /im spot, the sample temperature 
rise is estimated to be only 1 K. The spatial resolution 
needed for micro-Raman experiments is also facilitated 
by the fact that the resonance condition occurs in the 
UV (~ 3.8 eV) where the penetration depth of the ra- 
diation is small 10 nm). For a sample with a donor 
density of 10~^^ cm~'^ and a focal spot size ~ 1 /xm there 
are ~ 100 donors in the field of view. At this density the 
mean number of donor pairs with separations within 10 
nm is about 4. 

If resonance enhancement allows sufficiently large 
signal-to-noise ratio for micro-Raman scattering, an im- 
mediate question is what information of the phosphorus 
donors the Raman spectra carry, and how they are re- 



lated to donor positions and interdonor coupling. Here 
we focus on the Raman spectra of a pair of donors at 
and Rb (relative position R = R^i — R^), namely, when 
the |0} and |/) states of Eq. |(5J) are two-donor states. It 
is thus necessary to study the composition and energy 
spectra of the Raman-active donor-pair states, which is 
the theory we develop here. 

In the limit of distances jRj much larger than the ef- 
fective Bohr radii, consistent with the Si QC proposals 

1^ ^ , the two-particle singlet and triplet states within 
the IS manifold may be written in the Heitler-London 
(HL) approach: 

V:'^(ri,r2) ^ {u[\A{,Bi)±\BlA^)] 
t 

+v[\A'„B^,)±\B\,A'2)]}{\ n)TUT» (4) 

where 1^1,52) = V'i(ri— RA)V'i(r2-RB)- Here A and i? 
refer to the two nuclei locations, i and j refer to the single 
donor electronic states, 1 and 2 refer to the two electrons, 
and only one of the spin triplet is shown. Normalization 
and symmetry define the u and v coefficients. The two 
lowest eigenstates are singlet and triplet states obtained 
from the single donor ground state i = j = Ai. Their 
splitting J(R) defines the Heisenberg exchange coupling 
between the ground-state donor electrons [see inset (b) in 
Fig. Q , and is a direct measure of the degree of control 
achievable over the entangled electronic states given in 
Q. When the initial state |0) is a singlet or a triplet 
with i = j = Ai, the lowest Raman-active final states 
I/) arc obtained for i = Ai and j = E^'^\ where k = 1, 2 
refer to the two doublet states whose degeneracy has been 
lifted due to the axial perturbation of the donor pair. 

There are 2 initial and 8 final states for Raman scatter- 
ing from a donor pair in the IS' manifold. The 16 possi- 
bilities reduce to 4 Raman active combinations due to the 
absence of spin transitions (we neglect spin-orbit interac- 
tions) and the A ^ B symmetry of the HL state Q : Two 
of type {|0)„ !/(''))+} and two of type {\0)t,\&'>)t}, 
K = 1,2. The superscript ± indicates u — ztv in 

The Raman allowed lines are schematically represented 
in the inset (b) of Fig. ^ and their calculated values for 
a substitutional impurity pair positioned along the [100] 
direction are given in Fig.^ where triangles (squares) la- 
bel singlet (triplet) lines, and open (filled) symbols corre- 
spond to K = 1(2). The behavior of the calculated Raman 
shifts in Fig.^illustrates general trends obtained for any 
orientation of R, namely the overall convergence of all 
lines to the single-impurity value as |R| increases, super- 
imposed with an oscillatory behavior (more apparent for 
K = 1 for this particular direction) due to intervalley elec- 
tronic interference of Si conduction band minima 0, |^ . 
In general, Raman transition energies are anisotropic and 
very sensitive to the relative positioning of the donor pair. 

The propagation and polarization directions of the 
laser and Raman-scattered light, as well as the sys- 
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FIG. 1; Raman shifts of the four active transitions for a pair 
of substitutional donors placed along the [100] direction in Si, 
relative to the single impurity value Aq. Different symbols are 
defined in inset (b) and in the text. Insets: (a) Energy levels 
of a single P donor in Si leading to a Raman shift Aq = 105 
cm~^. (b) Schematic representation of the two-electron states 
of a P donor pair that participate in the Raman process. 

tern temperature, are controllable parameters in exper- 
iments. The usual notation to indicate the polarization 
scheme is P = kL(ffL,fjj^)kj^. Here we have consid- 
ered the following: Pi = Z{X,X)Z, P2 = Z{Y,Y)Z, 
P3 ^ Z([110], [110])Z and P4 = Z([110], [110])Z, where 
X, r, Z are the crystal axis [100], [010] and [DDI]. The 
relative Raman intensity of the different lines in a sin- 
gle spectrum is strongly dependent on the polarization 
components [see Eq.(I2l]. We also consider temperature 
effects in the populations of the initial states \0)s and 
|0)t, assumed to be in thermal equilibrium, i.e. 

nt/us = 3 eyip{- J/ ksT) . (5) 

In Fig. 13 we illustrate the combined polarization and 
temperature effects in the Raman spectrum for a pair of 
substitutional P donors along the [100] axis separated by 
~ 100 A(18 lattice parameters). Frames |2Ia) and|2Ib) 
are the data for low temperature (T = IK), while frames 
EJc) andEfd) are high temperature data (T = 50K). In 
addition, frames a) and|2Ic) refer to polarization Pi, 
i.e., along R. Here only two lines appear, because the 
K — 2 lines (solid symbols in Fig.QJ yield M = for Pi. 
The spectrum becomes more interesting for polarization 
P2, perpendicular to R, shown in frames |2b) and[5fd). 
Here all four lines contribute to the spectrum. At low 
temperatures [framesEJa) and[21[b)], practically all of the 
spectral weight is concentrated in the singlet lines, while 
for r = 50 K [frames |2Ic) and|2Id)] all lines contribute. 



FIG. 2: Polarization and temperature effects in the Raman 
spectrum of an impurity pair separated by ~ 100 A along 
[100]. Each vertical line represents a (5-function with weight 
(line intensity) proportional to its length. Frames (a) and (c) 
correspond to polarization Pi, parallel to R, and (b) and (d) 
to polarization P2, perpendicular to R. Temperatures are as 
indicated on the left. 



and the triplet lines are almost 3 times stronger than the 
corresponding (same k) singlet lines due to the triplet 
spin degeneracy. Polarizations P3 and P4 lead to Raman 
lines qualitatively similar to P2, but with reduced total 
intensity. We find that Raman intensities are stronger for 
polarization configurations such that fj^ and are either 
parallel or perpendicular to the relative position vector 
R (such as Pi and P2 in the present case of R along 
[100]), a feature that can be explored experimentally to 
determine the pair orientation. 

In Fig. 13 we present the individual line intensities for 
increasing temperatures. The intensity of same-K-pairs 
of singlet-triplet lines undergo a crossover close to the 
temperature J/fc^, a result easily obtained from Eq. jSJ- 
Thus J may be quantitatively estimated from the tem- 
perature dependence of the Raman line intensity. For 
J ^ ksTaiim the lowest experimentally achievable tem- 
perature, no effect should be observed by raising T, thus 
providing an upper bound for the value of J . 

We have also investigated the Raman signature of 
an ensemble of donor pairs, which is within the cur- 
rent experimental capability. We consider a sample with 
an ensemble of donor pairs of average relative position 
Ro — 100 Aalong [100]. For each pair, with a given R^, 
R^ varies over all possible fee lattice sites within a sphere 
of radius 10 Acentered at R^ + Rq. The resulting spectra 
are given in Figs. 0^a) and^fb). The Raman signature 
with polarization P2 shows an interesting T dependence: 
as T increases, the spectral weight shifts upward in fre- 
quency, so that a low-energy shoulder at P = IK becomes 
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FIG. 3: Relative intensities of the Raman lines of an impurity 
pair for which the exchange coupling is J = 0.47 meV. Note 
that there is a crossover between singlet and triplet pairs at 
a temperature close to J/fcs, allowing to extract quantitative 
information about J 
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FIG. 4: Raman signature of non-Poissonian distribution of 
donor pairs. In panel (a) and (b) the pair distribution has ax- 
ial symmetry, while in (c) and (d) the symmetry is spherical. 



a high-energy shoulder at T = 50K. Our results thus 
show that conventional Raman experiments may be use- 
ful in identifying whether non-Poissonian axial distribu- 
tions of donor pairs were actually achieved. Figs-Qfc) and 
0fd) give results for a different type of pair distribution, 
also of average relative position |Ro|, but with spherical 
symmetry: R^i varies over all possible fee lattice sites be- 
tween two spherical shells of radii |Ro ± 10 1 Acentered at 
Rs- Here the Raman signatures of the two polarizations 



are essentially identical, reflecting the sample's isotropic 
nature. 

In this Letter we have proposed and analyzed a pow- 
erful and novel optical technique for disentangling donor 
exchange coupling information in the Si QC architecture. 
Our proposed technique actually transcends the specific 
problem (i.e. Si QC architecture fabrication) that moti- 
vated us-the proposed micro- Raman spectroscopic tech- 
nique has the potential of becoming a versatile diagnos- 
tic tool in nanofabrication in general, since it can provide 
precise spatial information about relative positioning and 
entanglement of neighboring atoms in a host material, 
something that is very hard to do with any currently ex- 
isting tools in solid state physics. 
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